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f i ] On 1 2 March 2008, the Cassini spacecraft made a close 
encounter with the Saturnian moon Enceladus, passing 
within 52 km of the moon. The spacecraft trajectory' was 
intentionally -oriented in a southerly direction to create a 
dose alignment with the intense water-dominated plumes 
emitted from the south polar region. During the passage, the 
Cassini Radio and Plasma Wave System (RPWS) detected 
two distinct radio signatures: 1) Impulses associated with 
small water-ice dust grain impacts and 2) an upper hybrid 
(UH) resonance emission that both intensified and displayed 
a sharp frequency decrease in the near-vicinity of the moon. 
The frequency decrease of the UH emission is associated 
with an unexpectedly sharp decrease in electron density 
from ~90 el/cm’ to below 20 ei/cnr that occurs on a time 
scale of a minute near the closest encounter with the moon. 
In this work, we consider a number of scenarios to explain 
this sharp electron dropout, but surmise that electron 
absorption by ice grains is the most likely process. 
Citation: Farrell, W. M., W. S. Kurth. D. A. Garnett, R. E. 
Johnson, M. L. Kaiser, J.-E. Wahlund, and j. H. Waite Jr. (2009), 
Electron density dropout near Enceladus hi the context of water- 
vapor and water-ice, Geophys. Res. Lett.. 36, LI 0203, 
doi: 1 0. 1 029/2008GL037 108. 

1. Introduction 

[ 2 ] A major discovery of the Cassini mission was the 
observation of a substantial water geyser at the south pole of 
Saturn’s moon Enceladus [Porco et at., 2006; Waite et ai, 
2006; Hansen et ai., 2006 J. The emitted water 1) is believed 
to be the source of the large toroidal cloud of hydroxyl 
radicals with a peak density near 4 R s previously observed 
by HST [ Shemansky et at, 1993], 2) creates an extended 
neutral and plasma torus that drives processes in the inner 
magnetosphere [Richardson and Jurac , 2004; Jurac and 
Richardson, 2005, 2007], and 3) is the source of ice 
comprising the E-ring [Haff et ai., 1983: Horanyi et ai, 
2008]. The observation of both water and internal energetic 
activity at the moon places Enceladus on a growing list of 
astrobtologicai targets-of-interest that includes Titan, Mars, 
Europa, comets, and asteroids. 
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[3] Given this extraordinary find, the Cassini science 
operations team planned a set of close spacecraft passages 
by the moon to obtain further details of its structure and 
physical processes. One such encounter occurred on 12 
March 2008. Figure 1 shows the trajectory of the spacecraft 
m an Enceladus frame of reference where +x is oriented in 
the direction of the moon’s orbital velocity, +y is oriented 
towards Saturn, and +z is oriented out of the moon’s orbital 
plane. The water plume is illustrated in Figure 1 as the 
purple-shaded region approximately indicating the location 
of highest water density. The spacecraft moved primarily 
southward at .--14 km/sec, passing near and aligned with the 
water geysers. Given that Enceladus’ radius is ^250 km, the 
entire encounter (defined herein as ±5 radii about the moon) 
lasted for a little over 4 minutes, and the encounter with peak 
geyser activity itself transpired on a timescale of minutes. 

[ 4 ] The Cassini Ion and Neutral Mass Spectrometer 
(IN MS ) [ Waite et ai, 2004] found (on a previous encounter) 
that the geysers consisted of emissions similar in nature to 
cometary plumes ■ a dominance of water molecules but 
also the presence of methane, carbon monoxide, carbon 
dioxide, and both simple and complex organics [Waite et 
at., 2006]. During this 12 March 2008 encounter, the 
IN MS -measured water concentration was found to be 
<S O'/cm 3 before closest approach (CA) at 19:06:1 1 SCET, 
outside the plume. However, concentrations increased 
sharply thereafter peaking near 107cm 3 between 
'->19:06:30- 19:07:00 SCET. During this period, the space- 
craft was traveling aligned with the plume but also migrat- 
ing axtally-inward toward the gas center- intersecting its 
largest concentrations at approximately 300 km from the 
southerly-directed surface source. Passing southerly at larg- 
er radial distances, 1NMS continued to detect plume emis- 
sions at levels > 3 x lO’Vcm ’ as late at 600 seconds after 
CA. 

[ 5 ] As we demonstrate herein, the presence of large 
concentrations of water changes the very nature of the 
plasma-neutral interactions. Specifically, prior to CA, the 
plasma can be considered collisionless. In contrast, during 
the close encounter, the plasma electrons have many colli- 
sions with water vapor and ice such that its effect on the 
plasma inflow cannot be immediately neglected. The effect 
is to create a steady depletion of electrons with inward 
radial distance that is reminiscent of the exponential elec- 
tron decrease at the bottom-side of an ionosphere but now 
defined by water gas and ice chemistry. 

2. Observations 

[si Figure 2 shows a Cassini Radio and Plasma Wave 
Science (RPWS) radio spectrogram [Gumett et at., 2004J 
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Figure 1. The trajectory of the Cassini spacecraft past Enceladus, with (left) showing a ‘top’ view looking down from 
above the orbital plane and (right) showing meridional view looking opposite the direction of the orbital velocity vector, 
'fiie electron concentration drops below 70 cm ’ between the red markers. 


encompassing a 30 minute period centered about the Hnce- 
ladus encounter. Closest approach to the moon (52 km 
above the surface) is indicated in Figure 2. The measure- 
ments are from the wideband waveform system with a 3-dB 
bandwidth of ~7S kHz. The upper hybrid emission (f„ h ) and 
electron cyclotron frequency (f ce - derived from magnetom- 
eter results) are presented on Figure 2 for reference. 

[?] In this presentation, there are two dear and distinct 
signals: 1 ) A broadband signal (bandwidth of -“80 kHz) that 
commences near 18:55 SCET and extends to at least 
19:20 SCET that consist of impulsive spike-like events 
resulting from micron-sized dust grain impacts on the 
spacecraft system [Gurnett et at., 1983; Kurth et a!., 
2006; Wang et ai, 2006] and 2) and upper hybrid emission 


(labeled f nh in the spectrogram) that displays a. clear and 
distinct frequency downshift centered about CA, Unfortu- 
nately, the intense broadband dust impact signatures obscure 
the upper hybrid emission between 19:06 to 19:07:30 SCET 
when the impact rate is highest, making it difficult to track 
the band throughout the entire interval. However, a period 
of emission frequency down-shift is observed between 
19:05 19:06 SCET and up-shift to near pre-encounter 
levels between 19:07:30-19:08 SCET. 

[s] Regarding (1), water-ice dust grains are incident on 
the spacecraft and antenna, creating impact ionizations that 
are detected as a bipolar voltage pulse on time scales of a 
few milliseconds [Gurnett et at., 1983], When Fourier- 
transformed into a frequencv-versus-time spectrogram they 
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Figure 2, Cassini RPWS wideband waveform spectrograms showing both the broadband impulsive dust grain impacts 
and upper hybrid emission (f 11h ). The electron cyclotron frequency is also indicated for reference. Closest approach to the 
moon is indicated by CA. 
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Figure 3. The electron density (in cm 5 and m~ 3 ) as derived from the upper hybrid resonance emission. 


form a broadband impulsive event extending to the limit of 
receiver bandwidth. A clearly-presented derivation of the 
dust charge transfer transform is provided in the appendix of 
Wang et at. [2006], During the Enceladus encounter, dust 
grain impact rates progressively increased from 18:56 
SCFT, peaking on the order of 1000/sec between 19:04 
19:10 SCBT. impact rates progressively decreased thereaf- 
ter. During peak activity near CA, the bandwidths were 
observed to increase to large extents. The impact rate 
appears to have an abrupt decrease near 19:1.4 SCET, but 
still some level of impact detection persists intermittently to 
beyond 19:20 SCET. 

fy] Regarding (2), an upper hybrid emission is typi- 
cally detected throughout the entire inner magnetosphere 
[i Gurnett et at, 2005: Persoon et a!., 2005. 2006], It allows 
a derivation of the electron plasma frequency (since f pe = 
(f'uh fL)‘ /2 ) and thus electron density at a planet-centered 
radial distance of 2-10 Saturn radii. Well away from the 
moon, the upper hybrid emission (i.e., the electron density) 
varies in a progressive fashion, with variations on scale 
sizes of a fraction of a Saturn radii [see Persoon etal, 2006, 
Figure 5j. However, Figure 2 indicates that very near 
Enceladus, there is a significant reduction in i u! , indicating 
the presence of a clear and distinct “bite-out’ in electron 
density that was not anticipated - occurring on spatial 
scales of a few Enceladus radii, in fact due to pick-up 
processes [Tokar et at, 2006; Pontius and Hill, 2006] the 
plasma flow velocity near Enceladus has been measured 
to slow substantially from ~26 km/sec to <10 km/sec 
[Wahlumi et a!., 2005, 2008]. The plasma densities would 
thus be expected to increase by a factor of 2-3 in the near- 
vicinity of the moon in order to conserve flux. Clearly, such 
an increase is not observed for electrons and an unexpected, 
sharp decrease was observed instead. 

[To] figure 3 shows the derived electron density from the 
f uh emission which displays the dear and. distinct dropout in 
electron density. Note that between 19:05 and 19:06 SCET, 


the electron density decrease with radial distance from 
Enceladus, dnydr, is on the order of —80 el/rrf. This large 
negative gradient can be considered a primary' marker that 
defines the values of the source and loss terms in the 
electron continuity equation. By the lime Cassini exits the 
plume, the electron density has recovered to approximately 
its pre-encounter level. 

3. Interpretation 

[u] The electron continuity equation can be written as 
dajdt + V • (n e v,) = S - L where v e is the electron flow 
speed past the moon, S are electron sources (primarily 
photo-ionization) and L are electron loss processes that 
include electron recombination, electron dissociative attach- 
ment and a new term defining dust absorption. 

[ iz] Assuming a time-stationary situation (constant gas 
emission and steady spatial gradient) then dnjdt -■ 0. One 
can write an approximate 1-D electron continuity equation 
along the Cassini trajectory past the moon as: 

<9(n V v t )/ds ~ v e dttt/d& 4- n v (h e /(h 

= m»h 20 - M: - kdn H 20iL - %*««« 0) 

where the left hand side describes the divergence of the 
electron flux upon approach to the moon, and the right hand 
side accounts for sources and losses with a water photo- 
ionization source described by v «h 20 > electron recombina- 
tion losses with water ions (of assumed comparable density 
as electrons) described as k f nj, electron dissociative 
attachment losses described by k £ . n H 20 W, and electron 
losses via dust absorption defined as t/ d v c n e [Jackson et at., 
2008]. 

[is] We will now apply this equation during the period 
before CA between 19:05 19:06 SCET, during the RPWS- 
ofrsen ed very' steep electron density gradient of dnJds ~ 
-- 80 el/m 4 . The electron density, n e , was ~8 x 107m 3 at 
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the start of the drop-out period near 19:05 SCET and <2 x 
10 '7m 3 at the end of the period near 19:06 SCET (she hand 
is difficult to identify' between 19:06-19:07:30 SCET). We 
note that during this period, n H2 o was not at. its peak 
levels, but ranges from 0.2 1.0 x !u n: . 

[s4] Regarding the flux divergence during this time, both 
electron density and velocity are undergoing changes upon 
approach to the moon. Specifically, due to pickup/mass 
loading effects, the plasma flow speed was previously 
observed to progressively slow from corotationai values of 
26 km/sec at ~ 14 R^. to below 10 km/sec at 9 R e [Tokar et 
al., 2006; Wahlund et a/., 2005; Pontius and Hill , 2006; 
Wahlund et al., 2008}. Close to the moon, the flows should 
slow substantially (in theory, approaching zero). The previ- 
ously-observed plasma How measurements suggest a veloc- 
ity gradient that decreases with approach to the moon as 
Ovjds — 0.007/s. If there were no other sources or losses, 
the electron density gradient should become positive to keep 
the flux constant in the progressively slower flow. However, 
we instead observe a strongly decreasing electron density 
with inward distance, and conclude there must be a very 
substantial loss mechanism. 

[rsj Assuming that the decrease in flow speed described 
above continues info the plume, the flow speed would then 
become v e ~7 km/s near 4.3 R c . Using n e ~8x_ lO'Vnr. the 
divergence of the flux <7(n c v e )/<9s ~ -10 ft el/nr *s. The loss 
processes described on the right hand side of the continuity 
equation (1) must then account, for this rate on the left hand 
side of the equation. 

{ 16 ] The first term on the right-hand side of equation (!) 
represents the photoionization sources of electrons. The 
photoionization rate of H 2 0, t/, at Saturn is ~4 x 10 
sec [Richardson et al., 1 998] and for Cassini INiVlS mea- 
sured water levels at: ~10 H Vm 3 , we obtain a small source 
term, S, of +40 el/nr-s. 

[ 17] Electron dissociative recombination (H 2 CT + e -- H + 
OH) is considered the dominant recombination process 
[Richardson, 1992] and the process that appears to limit 
the molecular ion life time near the orbit of Enceladus 
[Sittler et al, 2008]. The recombination rate constant, k r , 
was calculated using ihe integral k r = (2e m) f oy(u) u 
f(u) du, where u is the energy in eV, and aJa) the cross- 
section of Mu! et al. [1983] varying approximately as ~ u ’ 2 . 
We assume a Maxwellian distribution of electrons with 
thermal energy u c , flu) * 2 ir~" 2 u e exp(~-u/u t ,}, such 
that f u""' f(u) du = 1. We find that k,~7 x 10 ,s nvVs 
assuming a 5 eV electron temperature. In a quasi-neutral 
situation, n H , 2 o+ ~ n e , and the losses from electron disso- 
ciative recombination, -k,ji£, are —45 e!/m 3 -s (as we 
discuss below, the plasma is most likely neutral to order 
of magnitude). As such, recombination appears to just offset 
photo-ionization sources during this period. 

[is] Losses via electron dissociative attachment can also 
be derived. Near 6.5 eV, an electron can dissociate water 
predominately into OH and H- {e + H 2 0 = OH + H-). The 
cross section for this process has a peak near 5 x 1{T~ nr 
and a width of Au 1 eV [ftakawa and Mason, 2002; 
Delory et al. , 2006}. The rate constant is thus k d ~ (2e/m) 0 : ’ 
cr^fUji) u d ffud) Au ~ 5 x 10 ' ' m7s, with u<j = 6.5 eV and 
u* = 5 eV. The loss rate for this process is thus — k,t n H20 n c ~ 
—40 el/m’-s, a relatively small value as well. 


[to] Herein then lies our dilemma: If we consider the 
usual processes associated with water molecules - photo- 
ionization, recombination, asid dissociation - along the 
Cassini transit, we find that none of these molecular source 
and loss processes at ~ a few 10’s of el/m 3 -s can account for 
the intense electron flux divergence at — 10 e el/ro-s. These 
water molecules and ion processes alone cannot account for 
the severe loss of electrons observed before CA. 

[in] However, as indicated in Figure 2, there is clearly the 
presence of water-ice (dust) in the environment and this also 
can alter the plasma character, in equation (1), the effect of 
this dust absorption is quantified via the —i] ( > v e n c term, 
where tfe = o d n d is the inverse of the dust/electron collision 
mean free path [Jackson et at., 2008]. The quantity n d is the 
dust density and < 7 d is the dust/electron collision cross- 
section, rr ajff, where a eft - is an ‘effective capture radius' 
about the dust grain. 

[ 21 ] We now consider the constraints on dust particulate 
absorption of electrons that are required to be consistent 
with the divergence of the flux (left hand side of equation 
(1)) - in essence reversing the argument to examine self- 
consistency, If we assume that the observed flux divergence 
is created solely and completely by dust, absorption, and 
define electron flux as g = n e v, we find that dg/ds ~ - 1 Cf el/ 
nr’-s, = - ija g. For n e = 8 x 107m 3 and v., = 7 km/sec, we 
find that r/d = g"’ dg/ds ~ 10 "/m. Given that tj d is also 
equal to cr d n d , we now want to get an independent 
estimates of dust density, n d . to determine if the self- 
consistent dust cross section is in any way realistic for this 
value of t/<j. 

[ 21 ] As indicated in Figure 2, during the Cassini Encela- 
dus flyby, dust (water-ice) was directly detected via impact 
ionization in an extended region about the body (±30 RJ 
symmetric about CA. Based on counting RPWS wideband 
waveform bipolar impact signatures, it was found impacts 
on the approximate order of 1 000 per second were detected 
between 1 9:04 19:10 SCET. So many impacts occurred 
during this period, that the automated counter run in the 
RPWS ground data processing system saturated. The inci- 
dent dust flux (#/second) being sensed by RPWS is F = n d v 
A > 1000/second. For a spacecraft speed of ~14 km/sec and 
effective spacecraft collecting area of ~0.4 nr [Kurth et al, 
2006; Wang et al. , 2006], this count rate translates to a dust 
density' exceeding ny > 0.2/m 3 . 

[ 23 ] Given the constraint on r; d and an independent 
impact-based estimate of n d , we find tin effective cross- 
section for the electron absorbers to be a ti -- try. n t ! < 5 x 
10 0 m 2 which corresponds to a predicted effective radius 
of less than about 1 millimeter consistent with electron 
absorption by small particles. Thus, the observed changes in 
electron flux can be considered consistent with the pres- 
ence of large concentrations (> 0.2/m 3 ) of small absorbers 
(<I mm). We note that these are bounding values and there 
may be larger concentrations of substantially smaller sub- 
micron sized grains that would go undetected via impact 
ionization with the RPWS antenna system [Kurth et al.. 
2006], The formalism also applies to water clusters (groups 
of molecules) that may be large in number and also 
absorbing electrons (although the process may be more 
similar to dissociative attachment at the cluster level), 

[ 24 ] In the above argument, the detailed electrostatic 
processes are effectively Tumped’ into the absorption cross 
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section, t7 d , which provides limited physical insight. We 
now consider art independent dust-plasma electrostatic 
model by Goertz [1989] that describes the micro-physical 
interaction between the plasma and dust. The mode! indi- 
cates that the electrostatic nature of the dusty-plasma is 
quantified by a P-lactor defined as (a u e ) (n d /n e ), where a is 
the grain radius and u e is the electron plasma temperature in 
eV [see Goertz, 1989, Figure 5, Table 1], For P > 10 12 m- 
eV, the larger concentration of grains creates overlapping 
Debye sheaths that leaves the inter-grain gas medium with 
an overall negative plasma potential [see Goetz , 1989, 
Figure 4). Applying the limits provided by the continuity 
equation (i.e., a e f»- ~ 10 " 3 m, n« 0.2/m 3 , u c ~ 5 eV, and n e 
~ 8 x 10' m 3 ) we find that P ~ 10 11 m-eVand conclude 
that the dusty-plasma should indeed be acting in a collective 
way. From Goertz [1989, Figure 5] at P 10 n m-eV, the 
plasma potential drop is predicted to be <!> ~ —0.7 kT e /e or a 
few volts negative, and should result in (at least) a factor of 
2 decrease in electron density. In fact, the RPWS Langmuir 
Probe detected a clear and distinct drop in plasma potential 
of a few' volts exactly coincident with the near-Enceladus 
electron dropout. 

4. Conclusions 

[25] We considered a number of possible loss processes 
to explain the near-Enceladus electron drop-out including 
water-ion recombination and water molecule dissociation, 
but found that these could not generate the sharp and 
distinct gradient. However, if we consider absorption of 
electrons via water-ice of sub-millimeter scale size, we find 
that indeed the observed water-ice concentrations are iarge 
enough to account for the electron depletion and associated 
plasma potential drop both measured by Cassini. 
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